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ABSTRACT 

Oil exploration in water deeper than 1000 feet has 
motivated the development of drilling platforms with 
dynamic responses superior to those of the conventional 
jacket structures. One of the most promising of these 
is the tension leg platform (TLP) , a semi-submersible 
platform held in place by tension members connected to 
the ocean bottom. 

This thesis evaluates the wave force transfer function 
for the TLP in heave, pitch, and roll motions. This is 
accomplished by assuming the transfer function for a 
single leg and then extending that result to include the 
effect of phase differences between forces on multiple 
legs, including the effect of wave spreading. 
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Cosine squared directional wave spreading is discussed 
and then dismissed in favor of a Gaussian distribution of 
wave direction. Using the Gaussian distribution it is 
shown that when the wave spreading has a standard deviation 
greater than 0.5 radian the exciting force transfer 
functions in heave, pitch, and roll converge rapidly to 
that computed for seas uniformly distributed over all 
incidence angles. 
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I. 



Introduction 



Present day platforms, in water depths from 10 to 1000 
feet, are primarily of welded steel jacket construction. As 
the present limit of about 1000 feet was approached many 
problems in fabrication and installation were encountered, 
indicating the need for a shift to a new technology. More 
significantly, with inceasing water depth the natural per- 
iods of jacket platforms increase from about one second in 
100 feet of water to approximately four seconds in 1000 feet 
of water. There is significant wave energy with a period of 
four seconds and dynamic amplification becomes a problem. 

The expected life of the structure is shortened due to high 
cycle - low stress fatigue damage. 

One alternative structure is a tension leg platform 
(TLP) . A TLP (see Figure 1-1) is a bouyant structure held 
in position by tension members secured to the ocean bottom. 
Because the platform is not rigidly connected to the bottom 
its natural modes of vibration are significantly different 
from a jacket structure of similar height. A TLP in a water 
depth on the order of 1000 feet will have natural periods 
in heave, pitch, and roll of about two seconds, well out of 
the region of high wave energies. 
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Since the TLP has only recently been accepted as a 
potential alternative for deep water drilling there are few 
studies of its dynamic response to wave forces. Such preli- 
minary studies are necessary to provide information for the 
design of the operating platforms. Response prediction re- 
quires specification of the modal exciting forces, for example 
in pitch, roll, and heave. For low everyday sea states the 
exciting forces and structure responses can be considered to 
behave in a linear fashion. Linearity of excitation and re- 
sponse suggests that a frequency domain solution technique 
is appropriate; especially when one considers that ocean 
waves are a random process most often described by a wave 
amplitude spectrum. This thesis presents an estimate of the 
modal exciting force spectra in heave, pitch, and roll in a 
way which accounts for structure geometry and wave spreading. 



8 




Figure 1-1 Tension Leg Platform 
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II. Problem Statement 



The modal exciting forces on a tension leg platform (TLP) 
will be dependent upon the frequency and the angles of inci- 
dence of the incoming waves. For some angles of incidence 
(at a given frequency) the forces on the various legs will 
act together to produce a larger total force or moment. For 
other angles of incidence the phase differences will cause 
the forces or moments exerted on the various legs to cancel, 
tending to produce a smaller total force. 

In this report the total modal force of the platform in 
heave, pitch, and roll will be divided into two parts. The 
first will be the force exerted upon a single, isolated leg 
by waves of the given frequency. The second will be a geo- 
metric term to take into account the effect of the phase 
differences between forces on each of the legs. This second 
term will be specifically calculated for a square TLP noting 
that the method could be used for any shape. 

The calculations so far imply a unidirectional wave 
amplitude spectrum. Rarely is this observed in practice. 

A reasonable wave spectrum can be approximated as being 
centered about one bearing from the structure and 
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o 

diminishing to zero at bearings 90 either side of the 
central bearing. Two types of wave spreading functions 
will be discussed. The first is the cosine squared spread- 
ing law which, although widely used, allows little if any 
leeway in varying the spreading to fit actual measurements. 
The second is the Gaussian distribution which does allow 
continuous variation of spreading simply by changing the 
standard deviation. 

The modal wave force spectra upon the structure will 
be the integral over all incidence angles of the directional 
wave amplitude spectra, weighted appropriately by the plat- 
form heave force, roll or pitch moment transfer functions. 
The resulting spectra will be normalized in this thesis by 
the results obtained when a uniformly spread wave spectrum 
is used. This is done to demonstrate the effect of wave 
spreading on the heave, pitch, and roll exciting force 
spectra. 
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Ill . Formulation of Equations 



A. Platform transfer function 

1. Single leg transfer function . The single leg wave 

amplitude to vertical force transfer function, is a 

s 

function only of wave frequency. Because the leg has circu- 
lar cross sections there is no dependence on 3 , the angle of 

incidence. This thesis assumes that r s (<o) is available, 
either experimentally or through another approximation method. 

2. Square platform transfer function . Consider plane 
progressive waves of frequency oj coming from an angle 3 off the 
platform axis. See Figure III-l. Each leg is thus subject to 
waves of the same frequency but different phase. For deep 
water waves the wavenumber is k = 2 jt _ = u 2 • Thus the wave- 

X g 

length is 




The phase shift in one wavelength is 2 tt . The phase 
shift, <{k, due to a distance, 1 ^, along the axis of prop- 
agation is thus 



<t>i 



2tt1^ 

X 




Consider the center of the platform as the point of zero 
phase shift. Let d be the leg spacing and a = d//2 be the 
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Figure III-l Regular Waves Incident on TLP 



distance from the center of the platform to the center of any 
leg. 

The distance causing a phase shift for leg 1 is 

1 = a cos(- - 6 ) . For leg 2 it is 1 0 = -a sin (7 - 8) . 

14 *■ 4 

Similarly for legs 3 and 4 the lengths are 1 - = a sin(- - 6) 

J 4 

and 1 ^ = -a cos(^ - 8) . From Equation III-l the phase 



shifts are: 



4> = - a cos (7 - 8) 

i g 4 



Note however that: 



4> a sin (7 - 8) 

2 g 4 



<0 . ,1T 

(j> s» - a sin ( - - 8) 
3 g 4 



2 

* =-- a cos (7 - 8) 

- g 4 



111-2 



cos(- - 8) = cos(-)cos(8) + sin(?)sin(8) 
4 4 4 

= (/2/2) (cos (8) + sin(8)l 



sin(-J - 8) =* sin(5)cos(8) - cos(-)sin(B) 
4 4 4 

= (/2/2) {cos (8) - sin(B) ) 



Letting T = cu 2 d/g and recalling a = (/ 2 / 2 )d then 
Equations III -2 become: 
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4> = - (cos(3) + sin($) } 

i 2 



T 



<j> =-- (cos(S) - sin(3) } 



2 2 



4> = ^ (cos(B) - sin(g)} 



s 2 



<J> (cos(B) + sin(S) } 



III-3 



3. Heave . The transfer function for heave for the 
total structure, is the sum of the contributions of 

each leg. Taking phase into consideration each leg contri- 
butes T s exp{i (cot+4>^) } . Thus the quantity of interest, the 
square of the magnitude of r^ t , is 



summing real and imaginary parts, and taking the magnitude 
squared of the exponentials this becomes 



Note however from Equation II 1-3 that <K = -4>i and 
$ 2 = -<t> 3 . Thus we have 




Incorporating the identity e 1 ^ 1 = cos(4>) + i sin(4>) 



1 r ht l 2 = 1 r s I 2 ({cos(<j» l )+cos(<t> 2 )+cos(<j) 3 )+cos(«J> lt ) } 2 
+ { sin (4>i ) +sin(4> 2 ) +sin(4> 3 ) +sin(4>„ ) } 2 ) 



III-4 
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cos(<t> : ) = cos( 4 > J sin(<J> x ) = -sin(<J> 4 ) 

cos(<j> 2 ) = cos(<J> 3 ) sin(<|» 2 ) = -sin(<{>j) 



and Equation III-4 becomes 



|r | 2 = |r s | 2 {2 cos($ x ) + 2 cos(<J> )} 2 



ht 



r ht | 2 = 4 | r | 2 (cos(«j> x ) + cos(4> 3 )} 2 . 



III-5 



By replacing <j> x and <j > 3 with their values from Equation 
III-3 and setting R^ = | r^ t l 2 /| T s | 2 this becomes 



JJI m 

= 4 (cos{-(cos $ +sin 8 )) + cos{|(cos $ -sin 6 )}) 2 

T T T T 

Rh = 4 lcos(- cos 8 )cos(^ sin 8 ) - sin(^ cos 8 )sin(^ sin 8 ) 

+ cos(s cos 8 )cos(^ sin 8 ) + sin(n cos 8 )sin(r sin 8)) 2 

22 z 2 



T T 

R. = 4 {2 cos(- cos 8 )cos(- sin 8 )} 
h -> 5 



III -6 



A form often more suitable for computer work is obtained 
by applying the identity cos 2 A = *5 + h cos 2A yielding 



r ht | 2 = 4 (1+ cos(T cos 8)) (1+ cos (T sin 8)) |r | 2 



III-7 



which can be implemented by recalling that T= <o 2 d/g. 
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4. Horizontal force assumptions . The roll and pitch 



moments experienced by a TLP will be dependent upon horizon- 
tal as well as vertical forces. In this report it has been 
assumed that the total moment producing effect of the hori- 
zontal force is small compared to the vertical forces. This 
is perhaps not always the case but is assumed in this discus- 
sion because, while the horizontal wave forces may be large, 
they are distributed along the legs, both above and below the 
center of rotation of the platform. The force on the leg 
near the surface is large but with a small moment arm. The 
forces below the center of rotation will be smaller but will 
extend further down the leg and thus have a large moment arm, 
producing a moment similar in magnitude but opposite in sign, 
and will tend to cancel the other. Thus the moment producing 
effects of the horizontal forces are assumed to be signifi- 
cantly smaller than those of the vertical forces. When the 
location of the center of rotation and the roll and pitch 
natural periods are known, as they would be for any particu- 
lar TLP design, the horizontal contribution to the roll and 
pitch moment could be calculated. It can be shown that the 
phase angles encountered when adding up the moments generated 
by the horizontal forces on each of the four legs would be 
the same as calculated previously for heave. Therefore the 
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results computed for heave will serve as a check on the qual- 
itative behavior of the role of horizontal forces in the 
generation of roll and pitch exciting forces. 

5. Pitch and roll . Pitch and roll for a square plat- 
form are exactly symmetric. Thus only roll will be discussed 
here with the understanding that pitch could be obtained by 
shifting the coordinates by 90 degrees. Taking the direction 
of positive roll moment to be as shown in Figure III-l the 
total transfer function for roll, squared, would be 

| T . | 2 = | r | 2 | (d/2)e i< *’ l + (d/2)e i '*’ 2 +(-d/2)e i<,>3 +(-d/2)e i<,>l *| 2 
rt s 

* I r g | 2 | 2 | e i< t’ 1 +e i( J , 2 -e i( J > 3 -e i< l ) '» | 2 

Let = I T rt | 2 / | T s | 2 . Incorporating the identity 
e 1 ^ = cos <{> + i sin <|> , summing the real and imaginary 
parts, and taking the magnitude squared this becomes 

d 2 

R = ^ {(cos <J> X + cos $ 2 - cos <J> 3 - cos <J>^) 2 

+ (sin 4>j+ sin <J> 2 - sin <t> 3 - sin <\) 2 } III-8 

However, from equations III-3 4> H = and 4>j = -<j> 2 . 

Thus cos 4>i» = cos , cos <f> 3 = cos 4>2 , sin <J> 4 = -sin <t> 3 , 
and sin <J> 3 = -sin <J> 2 . Therefore the cosine term in 
Equation III-8 reduces to zero and the sine term becomes 
(2 sin + 2 sin <J> 2 ) 2 , leaving 
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R = d 2 ( sin <$>!+ sin <f> 2 ) 2 . 



Incorporating the values from Equations III-3 this becomes 

T T 

R r = d 2 (sin{^(cos 6 +sin 8 ) }+sin{-- (cos B -sin B )}) 2 

T T T T 

R_ » d z (sin(-cos 8)cos(-sin B) + cos (-cos B)sin(-sin B) 
r 2 2 2 2 

+ sin(--cos B)cos(-sin B) + cos(-^cos B)sin(-sin 3)} 2 
2 2 2 2 

T T T t 

R = d 2 {sin(-cos 3)cos(rsin B) + cos (^cos 8)sin(Asin B) 
r 2 2 2 2 

mm T T 9 

- sin(-cos B)cos(-sin B) + cosUcos 8)sin(-sin B)) 2 
2 2 ^ 2 

R = 4 d 2 cos 2 (-cos B)sin 2 (^sin B) III-9 

r 2 2 

By applying the identities cos 2 A = h + H cos 2A and 
sin 2 A = H - H cos 2A this yields a form more suitable for 
computer work. 

| T rt | 2 = d 2 (1 + cos (T cos S ) ) f 1 - cos(T sin 8) ) 1 r s | 2 

III-10 

As mentioned previously pitch is simply rotated by 90 
degrees. Thus cos B becomes sin 8 and sin B becomes -cos 8 . 
The transfer function for pitch for a square platform is thus 

|r | 2 = d 2 (1 - cos (T cos B) Hi + cos(T sin B) ) | r | 2 . 
p t s 

III-ll 
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B. Wave spectrum . 



The wave spectrum incident on the platform varies with 
both frequency and direction and is denoted S (co, 8 ,B 0 ) . 

So is the center of an arc of 180° from which the waves 
are coming. 

For simplicity of calculations it is often assumed, as 
it shall be in this report, that the total wave spectrum 
can undergo separation of variables. Thus we can express 
the spectrum as 

S (w, 8 > 0 . ) = S (co) F(8 / 8. ) 

non o 

where the entire frequency dependence is contained in 
S n (co) and the entire angular dependence is in F(B,8 0 ). 

This section deals with approximations for the angular 
dependence , F ( 8 , 8 0 ) . 

1. Cosine squared wave spectrum . An often used method 
for describing the angular spreading of waves is the cosine 
squared wave spectrum (see Figure III-2) . It has the 
advantage that it vanishes for values of 8 equal to plus and 
minus tt/ 2. The formula for the cosine squared spectrum is 

— ^ p(S, 8 0 ) Sjj(to) $ -it/ 2 < 8-S 0 < tt/2 

where 8 0 is the center of the wave spectrum and C is a constant 
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used to normalize to unity the area between -tt/2 and tt/2. 
It is derived as follows: 

tt/2 

1 = / C cos 2 (6) dS 
-tt/2 

1 " /2 , 

- = 2 / cos 2 (8) d8 = 2 (7) 

C 4 



Thus the equation for the cosine squared spectrum is 
S (<o,8,8 0 ) = - S (w) cos 2 ( 8-8 0 ) ; ~ < (8-8„> < ? 

n M ir n 2-2 

III -12 

This spectrum suffers from one major disadvantage. 

It cannot allow for variations in the directionality of the 
waves. This problem can be eased somewhat by considering 
cosine to the fourth, to the sixth, and so forth but this 
results in only step changes in the directionality. Methods 
have been designed to provide continuous changes but they are 
somewhat cumbersome and comparison is not apparent. 

2 . Gaussian wave spectrum . The method used in this re- 
port to achieve continuous changes in directionality is to 
consider the waves to be spread about the central angle 8 # in 
a Gaussian fashion with standard deviation o. The reader not 
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-ir/2 



<B-M 



*/2 



Figure III-2 Cosine Squared Spreading 




Figure III-3 Gaussian Spreading 
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familiar with Gaussian (or normal) distributions may wish to 
refer to any good probability text for an explanation more 
detailed than the following. 

The probability distribution of a Gaussian process with 
mean Bo and standard deviation a is described by 

p(B,8 0 ) = — - — exp{- (8-B 0 ) z /(2a z ) } 

/7W a 

A spectrum of waves defined by a Gaussian process would 
thus be 

S n U,8,B 0 ) = C p(8, B 0 ) ; -ir/2<8-8 0 <ir/2 III-13 



where again C is a constant to normalize to unity the area 
between -tt/2 and tt/2. It is derived as follows, referencing 
8 0 to zero: 

tt/2 

1 = / C p(8,8o) dB 

-tt/2 

. tt/2 

± - / p(3, Bo ) d8 

C -tt/2 

This integral is evaluated using the properties of the 
Gaussian distribution and yields 



C = 



2 *(")-! 



where $(x) = 



1 x 

=-^= / exp(-t 2 /2) dt 

/2ir 
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is known as the normal distribution function and is tabulated 



in probability texts and math tables. Thus C is a function 
of the standard deviation and shall hereafter be denoted C(o). 

Unlike the cosine squared spreading the Gaussian distri- 
bution does not vanish for 8 equal to plus and minus tt/2. 

This is why C(a) varies with a. It must account for the 
cropping of the "tails" of the distribution for 1 8 1 > tt/2. 

Some simplifying approximations can be made however. 

As can be seen from Figure III-3 when o reaches a value of 
2 or larger the angular spectrum approaches a constant 
value. That value is 1/tt (in order to achieve unity area) . 
Thus Equation III-13 becomes 

S n («,8,8 0 ) = «* S n (w) ; ~ < 8-Bo < \ for o>2. 

For o < 1 the area in the tails diminishes rapidly 
leaving almost all of the area between -tt/2 and tt/2. As 
can be seen in Figure III-4 C(o) rapidly approaches unity. 

For o=1.0 C=1.13. For o=0.9 C=1.09. So for o < 1 
Equation III -13 can be approximated by 

S (to, 8, 8c) = S (<o) {“ 7 == — > e xp{ ~ ( ^~ 2 o) --> • 
n n /2tt o 2a 

For values of o between 1 and 2 the entire exact formula 
must be used 
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Figure III-4 Correction Factor for Loss of Gaussian Tails 
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S (u,B,B 0 ) s n ^ ^ 2 $ (I_) -1^/2tt exp ^ 
n 2a 



-(S-B n ) 

2o 2 



-} 



III-14 



3. Comparison of cosine squared with Gaussian . Since 
many people are familiar with cosine squared spreading it 
is useful to discuss for what values of a the Gaussian would 
approximate the cosine squared. Although the two functions 
are quite distinct there are several possible methods to 
compare them. 

The normalized cosine squared distribution has a value 
of 2 /tt * 0.637 for & = B»* For Gaussian spreading this would 
be C(o) / (/2 tt a) at B = Bo • If a is assumed less than 1 
then C(a) = 1. Equating the two values yeilds a = 0.627. 

When B = B 0 + tt/4 the cosine squared distribution is 
half of its maximum value. This situation can be matched 
with a Gaussian yielding still another "equivalent" a for 
the cosine squared. 



C (to) 
/2tT a 



exp{ 



-(k/4) 2 
2a" 2 



} = (0.5) 



CUO 
/2ir a 



exp { — — > = (0.5) 
2a 2 



„ r ~* (ft/ 4 ) 2 (h) ,h 

~ { In (0 . 5 ) } 



a - 0.667 
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The two values obtained are approximately equal. Thus the 
Gaussian distribution with c=0.65 is very similar to the more 
familiar cosine squared distribution. A smaller a means a more 
unidirectional wave front and a larger o means more spreading. 



C. Platform exciting forces . 

Each differential element of the spectrum S n (co,B,3 0 ) dB 
induces mean square force or moment upon the platform through 
the appropriate platform modal force transfer function, 
magnitude squared, |r t (co,B)| 2 . Thus the total force or 
moment spectrum is given by 

TT 

H(co,3o) = / |r (co, B) | 2 s (w, B , B o ) dB III-15 

— TT 

It is important to note that this modal force spectrum 
is dependent upon B 0 , the center of the directional spectrum. 

1 . Square platform, heave force . The heave force 
spectrum for a square platform in normally distributed 
waves is derived by substituting Equations III-7 and III-14 

into Equation III-15 yielding: 

8 0 +tt/2 

H (co,B 0 ) = / | T | 2 4 { 1+cos (T cos B) ) ( 1+cos (T sin S) )* 

B 0 -*/2 5 



S (co) 
n 



$ ( TT ) -1 
2c 



}{ 



/2tt a 



} 



exp<^£*i 
2 c 



-} dB 
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C(a) 

B h <».B,) - 4 l r s l S „ to) -jWl 



X 



S 0 + ‘ n '/2 in o \ 2 

/ {l+cos(T cos 8) } { 1+cos (T sin 0) } exp{^^4-} d8 III-16 

0 0 -i it/2 

Setting Iw(<o,8 ) equal to the integral in Equation 
III-16 and expanding the integrand yields 

0 o +tt/2 

1 (co,8 0 ) = / {1 + cos(T cos 0) + cos(T sin 8) + 

h B 0 ^/2 

cos(T cos 8)cos(T sin 8)) exp { - — ) d8 III-17 

where H.(a>, 8 0 ) =4 |r | 2 S n (o>) £££>. I. (<o,0 o > . 
n s /2ir o n 

A closed form solution of Equation III-17 does not 
exist. Thus a numerical integration was necessary 
( see Section IV ) . 

2. Square platform, roll moment . Similarly the roll 
moment spectrum for a square platform in normally distributed 
waves is derived by substituting Equations III-10 and III-14 
into Equation III-15 yielding: 
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B 0 +k/2 

H (oj , B o ) = / d 2 {l+cos(T cos S)}(l-cos(T sin 3)1 |r | 
3.-V2 S 



S n U) {_ i_) ( 1 ) expt^zilLl) 

n 2 /2ir o 2 

2a 



d8 



H (a,B 0 > = d 2 |r I 2 S (co) 7 Si q - ) - X 
r s' n / 2 tt a 



So+tt/2 



7 - { 8_g } 2 

/ {l+cos(T cos 8)}{l-cos(T sin 8)} exp{ — 26 2 ^ II1- 
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8 0 -ir/2 



Setting 1 ^( 01 , B 0 ) equal to the integral in Equation III-18 
and expanding the integrand yields 

8 o+tf/2 

l r (w,8 0 ) - f (l+cos(T cos 8)-cos(T sin 8)- 
8 0 -ir/2 



cos(T cos 8)cos(T sin 8)) exp{ ~ ^ 0 - — } 



2a‘ 



dB 



III-19 



where H r (co,8 0 ) 



d 2 



|r s l 2 s n <“> 



cw 

/2tt a 



I r (a» f Bo) . 



Again a closed form solution does not exist and the same 
numerical integration technique was applied ( see Section IV) . 
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3. Square platform, pitch moment . As mentioned 



previously, for a square TLP the pitch moment is simply 
rotated 90 degrees from the roll moment. Thus 



4. Normalization . In order to compare the effect of 
the wave spreading for different configurations a method 
of normalization was developed. The normalization factor 
was taken to be the force or moment frequency spectrum 
resulting from waves of equal amplitude from all directions 
in one half of the plane. Thus S n (<o,£,£o) equals (1 /it) S n ( cu ) 
for { 3 o — ir/2> < £ < (£ 0 +tt/ 2) and zero elsewhere. The 
normalization frequency spectrum becomes 



Because of the symmetry of a square platform and 
the uniformity of this spreading it makes no difference 
where the center of the spectrum is. For convenience 
£ 0 was taken to be zero and thus the normalization spectrum 
becomes 



H p (co , £ 0 ) = d 2 





H q (oj , £ 0 ) 



Qo+tt/2 

/ | T <03,6) | 

£ 0 -*/2 



2 (1/ir) S n (oj) d£ . 




tt/2 



III-20 
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which is simply the mean square value of T times S (co) . 

t n 

For heave and roll they are as follows 



4 */ 2 

H 0h^ = ^ S n^ |r g | 2 ^ Cl+cos(T cos 0)}{l+cos(T sin 0)} d0 

III-21 



H 



Or 




(co) 



tt/2 

| r I 2 / (l+cos(T cos 

S -tt/2 



0)}{l-cos(T sin 0)} d0 

III-22 



The normalized spectrum for heave is thus 



HN (co , 0 0 ) 
n 



H h (o) > 0 o ) . 
H 0h (< “ ) 
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Expanding the numerator and denominator and cancelling 

the common terms S (co) and If I 2 this leaves 

n 's' 



HN ( co , 0 0 ) 
h 



Ilfo (co , 0 Q ) 

!3 h (co) 
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where 

ll h (w,0 o ) 



H^ (co , 0 O ) 

s (co) |rl 2 
n 's' 



13 

h 



(co) 



4 C (q) 
/2tF o 



0O+TT/2 

I (l+cos(T cos 0) ) x 

0 O -V2 

{l+cos(T sin 0)} exp{ a - — } d0 

2o 2 




s n ( "> 




T 




tt/2 

/ (l+cos(T cos 0)}{l+cos(T sin 0)} d0 
-tt/2 
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Similarly for roll 



HN r (co,So ) 



H r (co,B 9 ) Il r (o>,B 0 ) 

H n (u) = 13 (w) 

Or r 



where 



II (<o,0 o > 
r 



H r (to , 3 o ) 

s n <“> lr s | 1 



2 So+tt/2 

_d_ C (q) / { 1+cos (T cos 0) } x 

/2? o 0 o -tt/2 



(l-cos(T sin 0)} exp{ — } d0 



13 (co) 
r 



H 



Dx. 



S n ( o») 



(to) 




d 2 n/2 

/ {l+cos(T cos 0 ) } { 1— cos (T sin 0)} d0 . 
* -it/ 2 
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IV. Numerical Prediction of Exciting Forces 



A. Discussion . 

Because of the complicated forms of the integrands of 
Equations III-17 and III-19 no closed form solutions to 
the integrals were found to exist. Numerical methods were 
employed to determine the variation of resulting forces 
caused by the variation in input parameters. Appendices 
I and II provide program listings and an explanation for 
their use. 

B. Assumed values . 

In order to use numerical methods of calculation 
values must be assumed for the relevant parameters. In 
this thesis the initial values assumed come from two 
Amoco papers published in the "Journal of Petroleum 
Technology." (See Reference 1.) These papers discuss 
a proposed vertically moored platform (VMP) which is 
identical to the concept of a TLP. Amoco' s VMP is proposed 
to have leg spacings of 160 feet and is estimated to have 
heave, pitch, and roll natural periods of approximately two 
seconds. Because of the popularity of the cosine squared 
wave spectrum a comparable standard deviation of 0.65 
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radians was chosen. 



The parameter of interest in all of the following 
analyses is HN, the ratio of the magnitude of the transfer 
function for Gaussian wave spreading to the magnitude of 
the transfer function for uniform wave spreading. 

When the Amoco TLP was analyzed in normally distributed 
waves with a 0.65 standard deviation the results were some- 
what surprising. The variation in HN from unity in heave 
was less than one percent for all values of wave spectrum 
incidence. Bo > (see Table IV-1) . A much greater variation 
of HN for roll was expected but HN still remained within 
ten percent of unity (see Table IV-2) . Clearly for this 
case an initial assumption of uniform wave spreading and 
the resulting simplifications in calculation of the force 
spectrum will provide very useful first pass information 
to the designer. 

However the Amoco proposal is by no means the only 
viable configuration for a TLP. One must consider the 
variation in HN as other parameters vary. 

C. Variation in leg spacing . 

The most obvious variation in a platform is the spacing 
between the legs. With the 160 feet of the Amoco platform 
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as a reference, leg spacings ranging from 100 to 250 feet 
were considered. Analysis results for heave are presented 
in Table IV-3 and for roll (and hence pitch) in Table IV-4. 

For all of these configurations HN for heave is still 
within one percent of unity. Roll suffers somewhat but 
the calculated HN stays well within 20% of one. For initial 
design work this accuracy is sufficient and the assumption 
of uniform spreading would be justifiable. 

D. Variation in period . 

Because the wave spectrum is continuous in frequency 
waves of all periods will be incident upon the structure. 
However, the period of paramount concern is the natural 
period of the structure. This varies from structure to 
structure and is somewhat controllable. Again taking the 
Amoco value, two seconds, as a reference, periods from one 
to three seconds were considered. The leg spacing was 
assumed to be 160 feet and the standard deviation 0.65 
radians. Analysis results for heave are tabulated in 
Table IV-5. The results for roll however vary enough to 
be visible on a polar plot of HN vs. Bo , the center of 
the wave spectrum, and are shown in Figure IV-6. When 
possible a plot has been preferred to a table because of 
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the greater amount of information compactly presented. 

For heave the variation in HN remains at about one per- 
cent reaching a high of two percent for a three second per- 
iod. For roll HN shows marked fluctuation but again stays 
within twenty percent of unity. While the natural period 
has significant implications with respect to the periods of 
significant wave energy it does not appear to significantly 
alter the modal force transfer functions as long as some 
directional spreading is present. 

E. Variation in standard diviation 

The final parameter to vary is the standard deviation. 
The standard deviation for uniform wave spreading is infinite 
while that for unidirectional waves is zero. However, as 
shown in Section III-B standard deviations above two radians 
approximate the condition of uniform wave spreading. There- 
fore almost all of the effect of varying standard deviation 
occurs between zero and two. One would thus expect a signi- 
ficant change in HN as standard deviation changes, particu- 
larly for small values. 

This expectation is correct. Standard deviation was 
varied from 0.125 to 2 radians for both heave and roll 
(D = 160 feet, Per = 2.0 seconds) . The results are plotted 
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in Figures IV- 7 and IV-8 for heave and roll respectively. 

Consider heave (Figure IV-7) . For standard deviations 
of 1.0 or 2.0 radians HN stays within one percent of unity. 
When a = 0.5 the variation becomes five percent. For 
a = 0.25 it exceeds twenty percent and when the standard 
deviation is 0.125 radians even the pattern of the curve 
changes. For roll (Figure IV-8) the results are similar, 
although somewhat more pronounced. 

However, if one can allow a 15% error in an initial de- 
sign approximation then the uniform wave spreading assumption 
is valid for wave directional spectra with standard deviations 
greater than approximately 0.5 radian. Note that cosine 
squared spreading is included in this category. 

F. Summary 

Numerous other cases were considered, simultaneously 
varying leg spacing, period, and standard deviation. In all 
cases where o >0.5 the resulting normalized amplitude- 
platform force transfer function had a magnitude within 
20% of unity. Thus for waves with any reasonable amount of 
spreading sufficient accuracy for initial design calculations 
can be obtained by using the mean square value of T(w,8). 

This result allows for a simpler initial calculation of 
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platform response to random waves in the manner of 
4 

Vandiver. For situations when a < 0.5 or when more 
accuracy is required Vandiver's method can still be 
employed using the computer programs of this thesis to 
determine HN (Vandiver's Cl). 
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SQUARE TENSION LEG PLATFORM, HEAVE EXCITATION 



WITH GAUSSIAN WAVE DISTRIBUTION 



LEG SPACING = 160 feet 

PERIOD = 2.0 seconds 

STANDARD DEVIATION = 0.65 radians 





Correction 


factor due to loss of 


tails = 1.016 




Factor due 


to uniform spreading 


(13) = 3.94 


So 




11 


HN h =H/I3 


0 




3.90 


0.990 


5 




3.91 


0.992 


10 




3.92 


0.994 


15 




3.93 


0.997 


20 




3.94 


1.000 


25 




3.95 


1.002 


30 




3.96 


1.004 


35 




3.96 


1.006 


40 




3.96 


1.006 


45 




3.97 


1.007 



Table IV-1 Normalized Platform Heave Exciting Force 
(D=160 feet, Per=2.0 sec, o=*0.65 rad.) 
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SQUARE TENSION LEG PLATFORM, ROLL EXCITATION 



WITH GAUSSIAN WAVE DISTRIBUTION 



LEG SPACING = 160 feet 

PERIOD = 2.0 seconds 

STANDARD DEVIATION =0.65 radians 

Correction factor due to loss of tails = 1.016 
Factor due to uniform spreading (13) = 0.235 E+05 



8, 


I 1 


HN r =H/I3 


0 


0.216 E+05 


0.919 


5 


0.216 E+05 


0.9 20 


10 


0.217 E+05 


0.922 


15 


0.218 E+05 


0.925 


20 


0.219 E+05 


0.933 


25 


0.222 E+05 


0.943 


30 


0.224 E+05 


0.954 


35 


0.227 E+05 


0.967 


40 


0.231 E+05 


0.981 


45 


0.234 E+05 


0.996 


50 


0.238 E+05 


1.011 


55 


0.241 E+05 


1.026 


60 


0.244 E+05 


1.040 


65 


0.248 E+05 


1.053 


70 


0.251 E+05 


1.068 


75 


0.254 E+05 


1.079 


80 


0.255 E+05 


1.086 


85 


0.256 E+05 


1.090 


90 


0.257 E+05 


1.092 



Table IV-2 Normalized Platform Roll Exciting Moment 
(D=160 feet, Per=2.0 sec, a=0.65 rad.) 
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SQUARE TENSION LEG PLATFORM, HEAVE EXCITATION 
WITH GAUSSIAN WAVE DISTRIBUTION 



PERIOD = 2.0 seconds 

STANDARD DEVIATION = 0.65 radians 

CORRECTION FACTOR DUE TO LOSS OF TAILS = 1.016 



LEG SPACING 


HN h 


HN h 


( feet) 


(min.) 


(max. ) 


100 


0.997 


1.003 


150 


0.998 


1.003 


200 


0.994 


1.006 


250 


0.995 


1.004 



Table IV-3 Normalized Platform Heave Exciting Force 
(Various leg spacings, Per=2.0 sec, o=0.65 rad.) 



SQUARE TENSION LEG PLATFORM, ROLL EXCITATION 



WITH GAUSSIAN WAVE DISTRIBUTION 



PERIOD = 2.0 seconds 

STANDARD DEVIATION * 0.65 radians 

CORRECTION FACTOR DUE TO LOSS OF TAILS = 1.016 



LEG SPACING 


HN r 


HN r 


(feet) 


(min. ) 


(max. ) 


100 


0.978 


1.027 


150 


0.933 


1.067 


200 


0.882 


1.114 


250 


0.846 


1.160 



Table IV-4 Normalized Platform Roll Exciting Moment 
(Various leg spacings, Per=2.0 sec, a=0.65 rad.) 
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SQUARE TENSION LEG PLATFORM, HEAVE EXCITATION 



WITH GAUSSIAN WAVE DISTRIBUTION 



LEG SPACING = 160 feet 

STANDARD DEVIATION = 0.65 radians 

CORRECTION FACTOR DUE TO LOSS OF TAILS = 1.016 



PERIOD 


HN h 


HN h 


(seconds) 


(min. ) 


(max.) 


1.0 


0.997 


1.001 


1.5 


0.996 


1.005 


2.0 


0.990 


1.007 


2.5 


0.994 


1.004 


3.0 


0.980 


1.024 



Table IV-5 Normalized Platform Heave Exciting Force 
(D=160, Various periods, c=0.65 rad.) 
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= 0 



HN = 



H 0r ( “> 



Figure IV-6 Normalized Platform Roll Exciting Moment 
(D=160, various periods , o=0.65 rad.) 
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Figure IV- 7 Normalized Platform Heave Exciting Force 
(D=160, Per = 2.0 sec, various a) 
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Figure IV-8 Normalized Platform Roll Exciting Moment 
(D=160 , Per =2.0 sec, various a) 
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V. Suggestions for Experimental Verification 



A. Single leg . 

Because the results of the computer program provide a 
transfer function between forces on a single leg to platform 
forces and moments, the single leg forces must be obtained. 
This can be done theoretically or experimentally. This 
report will deal with experimental determination. 

The single leg is assumed to be of circular cross sec- 
tion. Thus there is no angular dependence and testing can 
be performed in a ship model towing tank with a wave gener- 
ator. In the proposed test the single leg model will be held 
rigid with load cells measuring the forces as a result of 
incident waves. These forces will be normalized to forces 
due to scaled unit amplitude waves using a linear assumption 
for the force to wave amplitude relation. By performing 
this test with a variety of wave frequencies a plot of r s ( co ) 
can be obtained. From this the computer will predict 
platform heave force and roll moment transfer functions. 

B. Square TLP . 

To completely test the results of this study would re- 
quire the ability to generate a directional wave spectrum. 



47 



This ability is not available. However, the total transfer 
function, r (co,3), can be experimentally verified for heave, 
pitch, and roll. 

Once the results for a single leg are obtained a 
similar experiment with the total platform can be performed. 
The model should be held rigid with load cells measuring 
forces in all six degrees of freedom. The direction of 
wave incidence can be varied by turning the model. 

This experiment will generate a family of curves of 
r t (co , 3 ) , one curve for heave and one for roll for each angle 
of incidence 0. These curves can then be compared to those 
predicted using the measured forces for r s (co) and the heave 
and roll algorithms developed in this paper, Equations III-7 
and III-10. 
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VI. Conclusions 



A. Assumptions and approximations . 

1. Linearity . The forces computed in this paper are 
those exerted by waves on a motionless structure. The 
forces are assumed to be linearly related to wave amplitude. 
Furthermore the forces exerted on a motionless structure 
can only be useful in predicting dynamic response if the 
equations of motion are linearized. For a TLP the equations 
of motion will be linear for low to moderate sea states. 

This is due to the large dimensions of the major members and 
to the relatively small response amplitudes. Damping terms 
for the TLP may also be linearized because the damping 
ratios are expected to be small in heave, pitch, and roll. 

2. Effects of cross members . The model of a TLP used 
to derive the transfer functions had four legs acted upon by 
wave forces, and nothing else. Real TLP ' s must have cross 
members between the floatation legs to provide structural 
support. These would also be subject to wave forces. 

The diameter of a cross member will be much smaller than 
that of the leg. In the Amoco platform referred to in Refer- 
ence 1 the cross members have an outer diameter of 5.5 feet 
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while the floatation leg's diameter is 30 feet near the sur- 
face where the large wave forces are. This is a factor of 
about 5h difference. The displaced volume is proportional 
to the cross sectional areas which are almost a factor of 30 
different. In addition the support members are mostly near 
the bottom of the legs, away from the largest wave forces. 

The parts that are shallow are close to the legs and would 
react in phase with them. Thus the effect of the waves on 
the cross members can be assumed to be negligible with 
respect to the effect on the legs. 

3. Wave disturbances of the legs . Another aspect not 
considered in this model is the effect of one leg upon the 
waves incident upon the others. In the Amoco example the 
outer edges of the legs at the water surface are 130 feet 
apart which is only about four times their 30 foot diameters. 
There could be some shielding or reflecting effects present, 
causing the results of this study to be in error. As the 
ratio of leg spacing to leg diameter increases these errors 
will decrease. 

The primary purpose of this thesis was to evaluate the 
influence of wave spreading on the heave, pitch, and roll 
force spectra. The most important geometric effect is the 
relative phases of forces exerted on the legs. In cases 
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that the wave spreading is shown to be sufficient to average 
out the relative phase effects, then it can certainly be 
argued that interaction effects between legs would also 
be suppressed. 

B . Utility of results . 

1. Initial design approximation . For wave directional 
spreading with a standard deviation of greater than 0.5 
radians the heave force or roll or pitch moment imposed 
upon a TLP can be approximated within 20% to be that imposed 
by uniformly spread seas. By incorporating this into the 
method used by Vandiver^ an initial estimate of the response 
of a TLP can be obtained which would be applicable in 
preliminary design. 

2. More accurate approximation . If the standard 
deviation of the wave spectrum is less than 0.5 radian or 
if an answer with more than 20% accuracy is required then 
the computer program outlined in Appendices I and II can 
be used to determine HN (Vandiver's Cl). This can then 

be incorporated into the calculations of the modal response 
of the structure. 

3. Variation to fit other structures . Results for 
other square TLP's can easily be obtained simply by changing 
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the values assigned for the leg spacing and the period of 
interest. If the structure of interest is not square the 
same concept can be applied by changing the program defini- 
tions of the transfer functions. In this manner results 
can be obtained for any shape of TLP (square, triangular, 
rectangular, etc.) with any leg spacing and any natural 
period. 

4. W ave directional spectrum variation . The greatest 
advantage of modeling the wave directional spectrum as a 
Gaussian process is the continuous variation in spreading 
it provides. The standard deviation of the wave spectrum 
at a proposed location for a TLP could be measured under 
various conditions, noting that the periods of the signi- 
ficant wave energy would also change. This information, 
standard deviation and period, could then be input to this 
program to estimate platform forces and moments under those 
conditions. The final information of forces, moments, and 
responses would then be applied to sizing the tension 
members for adequate tensile and fatigue strength and to 
ensure the platform motions do not cause equipment malfunc- 
tions or personnel discomfort. 
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NOMENCLATURE 



a distance from center of TLP to center of any leg 

C constant to normalize the area of the wave 

directional spectrum 

d leg spacing of TLP 

g acceleration of gravity 

h as a subscript, denotes heave 

H x frequency spectrum of the total force or moment 
on the TLP due to wave spreading 

H Qx frequency spectrum of the total force or moment 

on the TLP due to uniformly distributed waves 

HN X normalized frequency spectrum of the total force 
or moment on the TLP 

I integral portion of H dependent on wave spreading 

X X 

Il x numerator of HN x , dependent on platform geometry 
and wave spreading 

I3 x denominator of HN x , dependent only on platform 
geometry 

k wavenumber 

1^ distance from center of TLP to projection of leg i 
on the direction of wave propagation 

p as a subscript, denotes pitch 

r as a subscript, denotes roll 

R ratio of magnitude squared of platform T to 

magnitude squared of single leg T 

s as a subscript, denotes single leg 

S n wave spectrum 
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t as a subscript, denotes total platform 

T non-dimensional quantity determined by leg spacing 
and wave frequency 

S angle of wave incidence measured from roll axis 

B 0 angle of incidence of center of wave spectrum 

T x transfer function from incident wave amplitude to 

force (or moment) exerted upon the structure 

5__ Dirac delta function 

mn 

X wavelength of incident waves 

a standard deviation of wave directional spectrum 

Z summation 

phase shift of the wave incident on leg i as compared 
to the wave at the center of the TLP 

$ normal distribution function 

w radian frequency of incident waves 
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Appendix I . Computer Program Listing 



A. Discussion . 

Because of the complicated forms of the integrands of 
Equations III-17 and III-19 no closed form solution to the 
integral was found to exist. Yet the trigonometric func- 
tions with trigonometric arguments caused such rapid oscil- 
lations of the integrand that most numerical methods failed 
to converge. The method found to be the most useful was 
Gauss-Legendre quadrature (see Appendix III) . A quadrature 
of fifth order was found to give sufficient numerical 
accuracy and is contained in the subroutine GAULEG. 

Because the wave directional spectrum was only non- 
zero between plus and minus tt/2 radians a correction factor 
was necessary to ensure its area was unity. This correction 
was accomplished by an eighth order polynomial approximation 
to the area under the normal curve. This approximation was 
developed by the author and is contained in the subroutine 
ERF. 

Because of peculiarities in the computer system 
available the platform transfer functions are entered as 
function subprograms HEAV and FHEAVE for heave, and ROLL 
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and FROLL for roll, the former of each pair being the 
transfer function with uniform wave spreading. 

The program is written in FORTRAN-IV language and has 
numerous comment statements to aid in its understanding. 

B . Program listing . 

A listing of the programs for heave and roll, with 
associated subprograms, follows. 
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C I Cl C ) t) MCI MCI MCI 



TENSION LED PLAT^OR* 

u?i vg 

EXTERNAL wr AV,FHEAV 

SET STAND? RO DEVIATION (SD),LSG SPACING ( D ) , AN D PERIOD (PER) 
r>=i so. 

S0:O.55 
P r R =2 . 

CALCULATE VARIANCE AND t 
200 VAR-5D**2 

300 T-1.226C3B*D/PFR**2 

W»I*t (5,5) 

5 r 0 ° ~ A T ( '1 *,ri~, •SQUARE TENSION led platform, REAVE EXCITATION*) 

writ: (5, ID) 

10 FORNAX (• * , XI 5 , * C A0S3 XA N WAVE DT STRI BU TI ON • //) 

V°IT" (5,15) D 

15 FORMAT (’ *, T 15, , LES 5? * CTND = • , m Ud , FS . 0 , • FE r T # ) 

WRITE (5,20) P rD 

20 r D c XAT (• ^T**, *PEFIDD=* ,T4D,^.1 , * SECOND^*) 

WRITE (6,25) 5 D 

25 *0**kT C * , ?1 c , • STAND A°0 D E VI A TIDN* • , T ?3 , F r . 2 , * RADIANS*//) 

CORRECTION' FACTOR FDR LOSS OF TAILF OF CAUSSIAN 
WIDTR = 1 .^D7n5/SD 
CALL ERF (WIDTH, PHI) 

C = 1 ./(2.*P U I-1 . ) 

P! Tr ( 5,30) C 

3} FOP M * T C *,115,* DIRECTION FACTOR DUE TO LOSS OF T AILS=* ,F6 .3 ) 
FACTOR DUE TO UNIFORM S D R?.ADI*7 

XL*-1. 570756 
X'J- 1.570796 
d ^ - d . 

CALL DAUL^C (XL,XU,H r AV,R,I p , n ^,VAP,T) 

PI 3 = P 

WRITE (5,35) R 1 3 

35 POP w 5 T (* *,116,* FACTOR DUE *"0 UNIFORM S^RE * DTNC ( 13 ) = * , El 1 . 4 ) 
V"I"E (5,40) I R 

40 FORMAT C * ,T2 r , * (ERPD° CONE 1110*= * , 13 , * ) * / ) 

WRI^E (5,45) 

45 FD D * A*" (* *,* BETA* ,123, *11* ,?4 5 , # KN = I1/I3*,T52,*SR50R?*/) 

V A ^ V DI or cnov OF CENTER O r WAVE S^ECTRU* 

?D0FD = '> . 

100. ED^SODFDM .7U5329E-2 
XL1*ZL+90 
V'M-VU + BC 

CALL DAULDC C X L 41 , X U 1 ,*H r * V , V ,IEP , PC , V AP , T) 

?7i -Y*C/ED 
H:I = RI1/PI3 

WRIT? (6,50) BODED, R T ^, p N,TEP 
50 FD OM A T (• ",F5.1 ,~17,S11.U,TU2,F10.4,T62,I3) 

P jDED=BDD"D+5. 

IF (PODED .LF. 45.) ~D TO 100 

STOP 

END 
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no on noon 



TENSION LEG ruATFORM 
ROLL 

EXTERNAL ROLLtFRCLL 

SET STANDARD DEVIATION < S D ) * L E 0 SPACING ( 0 ) * A NO PERIOD (PER) 

D= 1 6 0 • 

S D = G • 6 5 
PER=2. 

calculate variance AND T 
2 ( C VAP=$D**2 ‘ 

300 T=1 .22603fi*C/PER»*2 

WRITE (6*5) 

5 FORMAT ( *1**T1Q**SGU/RE TENSION LEG PLATFORM* ROLL EXCITATION*) 
WRITE (G* 10) 

10 FORMAT (* ♦ *T15* •GAUSSIAN WAVE DISTRIBUTION*//) 

WRITE (6*15) 0 

15 FORMAT (* * * T 1 5 ♦ *LEG SP AC I NO= * * T4 0 *F5 . 0 * * FEET*) 

WRITF (6*20) PER 

2 C FORMAT (* **T15**PERI0r = *tT*C*F5.1** SECONDS*) 

WRITE (6*25) SD 

25 FORMAT (* * *T15* ’STANDARD OE V IA T I ON= * * T2 9 *F6 . 2 * * RADIANS*//) 

CORRECTION FACTOR FOR LOSS OF TAILS CF SA'JSSIAN 
U I D TH = 1 .570796/SO 
CALL ERF (WIDTH*PHI) 

C-l • / ( 2» » PM I •! • ) 

WRITE ( 6 * 3 G ) C 

30 FORMAT (* * • T 1 5 * *C OP RE C T I ON FACTO* DUE TO LOSS OF T A I L S = * * FG . 3 > 
FACTOR DUE TO UNIFORM SPREADING 

XL=-1* 570796 
XU= 1.570796 

e*=o. 

call qauten (xl*xu*roll*r*ir* n O*VAR*T) 

RI2=R*D*+2 
WRITE (6*35) P 1 3 

35 FORMAT (* * * T 1 5 * *F A C TOR CUE TO UNIFORM SPREADING ( 13) = ** Ell. 4) 
V*ITE (6*40 IR 

40 FORMAT (* * * T20 * * ( r R PO R COND I TI 0N= * * I 3 /> 

WRITE (6*45) 

4^ FORMAT (• *** FETA * *T?3 ♦* 1 1 * *T45 * *HN=I 1/ I 3** T62* * ER ROR? •/) 

C 

C vary DIRECTION OF CENTER OF WAVE SFECTR'JM 
8 CDEG = 0 ♦ 

100 P0=90CEF*1.74532PE-2 

XL1=XL*30 
XU1=XU*B0 

call caulr^a (xll*xul*frolt*y*irr*DO*VAR*T) 

R I 1 = Y*C**2*C/SD 
HN=RI1/RI3 

WRITF (6*50) EWDEG*RI 1*FN*IER 
50 FORMAT (* **F5.1*T17*EU.4 f T42*F10.4.T62* 13) 

B0nEG=D0DEG^5. 

IF (80OEG .LE. 90.) GO TO ICO 

STOP 

END 
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SUa?OCJTTv r ( XT ,X 0 , r CT , T , IF* ,*0 , V? © , T) 

C THIS IS ' j?i ttu ORPEP SA r JSS-L?3EV0PE QUA D C A T T JP E 

1 = 0 

H = (VTI-XL)/i:. 

A1= C.23^2^ 

*2=0.47962? 7 
A3=C. 56999*9 

r 

X1 = f>. 9061793 
X2 = 0. 533469 3 

sux-o. 

c 

1 u=.5*^ 

r T2=.5*H 
A = X l+ fJ ? 

T 1=K1*P2 
Y2=X2*K? 

1 = 1 + 1 
sn*i=su* 

S0*=9. 

3 B 1 = * + T 1 
92= A- Y 1 
B3=A+*2 
P4=»-Y2 

SUX=SCTK*A1*( 7 :i(31 ^0^/AP # T)+^07(B2 ,90,7*9,11) + A. 2* 

1 (FC7( B3 ,90 ,T)+FCT( 94 , BO , V* * , T) )+A3*FC^( * ,R0 ,VA D ,T) 

C 

A = A+H 

IF ( XU-A )4,u ,3 

4 SU w =F?*SVi 

if rstn .F2- o.) so to * 

£RP=ABS((SU V -SU V 1) /SU M ) 

C 

IF fl.GT.IC) TO TO 6 
IF (FRR-.'M) 5, c ,1 

= , I,,., 

Y=S"* 

RETURN 
6 TFR = 1 
Y = SO" 

FFTHR’? 

EVD 
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; r ?^ 



'!'•! r u 



1 

2 



1 o 



7) T '* 

T —r 

t- “ v/ 

T r» r 

1 = 1 

T V 

±. 

PHI 

CO 

V = 1 

? :J I 

A r T v 
A "* v 
2, ? 7V 
A 'TV 

MJX 
A TT v 
A ’T v 

2 rT V 

no 

V = V 

?4I 
I ? 

7 r t 

r '-:c 



TT' C 'n 7 *’ f ’ ’7 7 7 f V r » ~ r ^ 

• N t % pf'T y vpv T ' y ■** t 7 *»j r 

^, r . „ ^•rnp P *A r ^ p *r a r ~TT o 1/ r 
^ ,7 oj r ” 2 »j v (o) 

a t r 3 »/ c 'a d ‘|rr? *t» t y t? 7 * t M P c 0 c v 

v ) 1,^2 



1 r- r *, 



’^0^7 



r Y 
= 1 



. L ' 



3,°> 



- 

f 1 
(2 
(3 
( u 
( = 

( 

( 7 

r 2 

1 r 

*Y 

= ?-TI fY*£T'f I 1 



1= 

)=-1.2^356—u 

) =- 6 . 1 ^ 127^-2 

) = - 1 • 2 0 7 j 4 r r - 2 
)= ^^c^^qor.o 

) = -0 .0 s — 3 

) = i,43? 3 Si5-3 
)=- 3 , 63 1 52 ? ,? -^ 
1 = 1 ,* 



- {- 

(T 

? TP 



El, 1 ) ? H T = 1 . -rcr 
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FUNCTION F nE A V < 3 t BO ♦ , T ) 

FHEAV=l.f9*76V*{ l.+C PS' T *C 0? (*>>>*< l C E • T » S I < _ ■ ) ) 
! E X < U - H 0 > * * 2 / ( 2 • * V A ” ) ) 

RETURN 

END 



function H E A V { L- * 2 0 ♦ V a " * T } 

HE A V = 1 * 2 7 3 2 A * t 1 . *■ C 0 S f T*CCS ( R) ) ) * ( 1 , * C 2 C ( T * S I . < 3 ' ) * 

RETURN 

END 



r U . i 2 T 1 3 N’ "P^EL ( r - , n 0 , V * p ,T> 

~ROLL=. 3-5F3U2* f 1 . + (E ) ) 1 * ( 1 .-''PS (T*S T '! (Pi ) ) * 

1 E v ^f-(°- t3 n**2/( P. *V'P) ) 

3 m 

F ” C 



F ij V I FT T C l) ?rL(",nr fV ! 0 f T) 

p 1LL= .3 1 f R 3 30 * ( 1 . +C n S( T* r PS f IJ ) ) ) *( 1 . - "’PE f T* 5 ( B ) ) ) 
° 2 T U R '! 

END 
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Appendix II. Directions For Use of Program 



A. Discussion . 

This program was intentionally designed to remain as 
much as possible in discrete sections, each performing a 
specific purpose. In this manner the program can be more 
easily understood by others and converted to use for many 
different sizes and types of structures. 

An error flag has been added to the subprogram GAULEG. 

If GAULEG is unable in a reasonable time to get consecutive 
iterations to agree within one percent it flags this fact 
to the user by indicating an error condition of one. If the 
error condition is zero then GAULAG has gotten the same an- 
swer (within one percent) for two iterations with the 
second iteration having twice as many increments as the first. 

B. Numerical parameters . 

To make it easier to understand the program flow the 
leg spacing (D) , standard deviation (SD) , and period (PER) 
are simply assigned values in FORTRAN statements. They could 
just as easily be read from data cards, input from an inter- 
active terminal, or varied in do-loops, whichever is most 
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convenient for the user on the available system. 

To adopt the program to a square TLP of interest one 
must simply replace the assignment statements for leg spacing 
and for the period of interest. If many periods are import- 
ant then a do- loop for period should be used returning to 
statement 300 to ensure the value of T is updated. If a 
do-loop is used to vary standard deviation it must return 
to statement 200 to revise the value of the varience. 

The center of the wave spectrum (3o) was varied in an 
implied do-loop with statement 100 because one is normally 
interested in the response to waves coming from all directions 
Because of the symmetry properties of a square TLP the angle 
must only be varied over 45 degrees for heave and 90 degrees 
for roll. 

C. Conversion to other types of platforms . 

The information concerning the platform geometry is 
contained in the platform transfer functions. In the computer 
program the definition of the functions in the function 
subprograms (ROLL, FROLL, HEAV, and FHEAV) entirely defines 
the structure. Note that a factor of d 2 was left in the main 
program rather than ROLL and FROLL since this factor will be 
present in moment transfer functions for any type of structure 
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Thus to change this program to accomodate any shape TLP 
one must simply perform derivations similar to those in sec- 
tion III-A and use the results to replace the function defin- 
itions in the function subprograms. 

D . Conversion to other wave spectra . 

This program was explicitly written for a normally dis- 
tributed wave directional spectrum because it is so easily 
varied. Simply by varying the standard deviation one can 
obtain as directional or uniform a spreading as desired. 
However, if another spectrum is preferred this program can 
be adapted. 

The most obvious place to change is in the definition 
of FROLL and FHEAV. There the new spectrum would be substi- 
tuted for the exponential term. In the main program the 
correction factor C/SD is applied after calling GAULEG using 
FROLL or FHEAV. This would be changed or eliminated. Since 
C, the factor for loss of tails, is no longer necessary the 
entire subroutine ERF could then be eliminated. New correc- 
tion factors for the new spectrum would then be introduced. 

It is also possible to normalize the transfer function 
using something other than uniform spreading. For this ROLL 
and HEAV only must be revised. If correction factors are 
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necessary they are best introduced in the main program after 
calling GAULEG for the integration. 
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Appendix III. Gauss-Legendre Quadrature 



A. Legendre Polynomials 

1. Definition . Define a set of polynomials P n (x) of 
degree n, where n is a non- negative integer, such that 



_{ P m (x) V x > dx = c n 5 mn 
and P n (l) = 



The first is an orthogonality condition and the second 
is necessary to fix c n . Together they define a unique set of 
polynomials known as the Legendre polynomials, the first few 
of which are 



P^x) = x 
P 2 (x) = % (3x 2 -1) 

P 3 (x) = h (5x*-3x) . 



In general they can be found by 



P (x) = — —-fj- (x 2 -l) n . 

n 2 n n! dx n 



Note that P (x) has n zeroes in the interval from minus one 
n 

to one and is even for n even and odd for n odd. 
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2. Polynomial approximation . Any function piecewise 



continuous in the interval from minus one to one can be 
expanded in terms of an infinite series of Legendre polynomials. 
It can be approximated by a finite series of Legendre poly- 
nomials. 



y (x) 



00 

* T. c 

m=o 



m p m* x * “ m I 0 c m p m* x) 



B. Integral approximation 

A simple means to evaluate the integral in the region 
minus one to one is as follows: 

1 1 00 

/ y(x) dx * / Z C-, P ra (x) dx 
-1 -i o 

Multiply both sides by 1 = P Q (x): 

/ y(x) dx ■ I c } P_ ( x ) P_(x) dx 
_ j 0 m - 1 m o 

1 

= c Q / P Q (x) P Q (x) dx 
} y(x) dx =* 2 c Q 

Thus the problem has now become one of finding c q , the 

coefficient of P (x) . One possible method will be shown 

o 

using the example where n=2. 

We desire to find coefficients and arguments such that 



68 



_} y (x) dx = 1 c 0 = A 0 y(x # ) + A x y(x t ) . 

For any two points. x 0 and Xj , we expect 

c . p , <\> + c , p , < x ,> + c 2 p 2 <\> - y<\> 

c . p . <x .> + °, p , tx ,» + c i W - y< x ,> 

Noting that P Q (x) = 1 and P t (x) = x and choosing x fl and 

x, to be the zeroes of P (x) this becomes 

c o + c x x « - y<V 

c o + c i x i = y (x ,> 

Where x # and Xj are ± (1//3) . By adding these two 

equations and noting that x Q = -x t we have 

2 C o = y(-l//3) + y(l//3) 

which is our desired answer. 

To get a larger degree of accuracy a larger n is required. 
This method for determining 2 C 0 remains effective. For the 
case of n = 5 the solution is (to seven significant digits)^ 
i 

(y (x) dx = (0.568 888 9) {y (0) } 

+ (0.478 628 7) {y(.538 469 3)+y(-.538 469 3)} 

+ (0.236 926 9) (y(.906 179 8)+y(-.906 179 8)} . 
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However, one is not always interested in the interval 
from minus one to one. A transform must then be used to 
transfer the region of interest to that interval. For addi- 
tional accuracy without the need for a larger order quadra- 
ture the total interval can be divided into subintervals. 
Each subinterval can be transformed seperately and the 
result summed. Thus each subinterval will be approximated 
by a fifth order polynomial. This is the method employed 
in the subroutine GAULEG. 
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